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TEMPERATURE HISTORIES IN CERAMIC -INSULATED HEAT-SINK NOZZLE 
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SUMMARY 

Temperature histories were calculated for a composite nozzle wall by 
a simplified numerical integration calculation procedure. These calcula- 
tions indicated that there is a unique ratio of insulation and metal heat- 
sink thickness that will minimize total wall thickness for a given operat- 
ing condition and required running time. The optimum insulation and metal 
thickness will vary throughout the nozzle as a result of the variation in 
heat-transfer rate. The use of low chamber pressure results in a signifi- 
cant increase in the maximum running time of a given weight nozzle. 

Experimentally measured wall temperatures were lower than those cal- 
culated. This was due in part to the assumption of one-dimensional or 
slab heat flow in the calculation procedure. 


INTRODUCTION 

One approach to providing exhaust nozzles suitable for long-burning- 
time, high-performance solid-propellant rockets is to insulate the nozzle 
with a high-temperature ceramic coating. Successful ceramic coatings 
have been the object of much research and development over the past 10 to 
15 years. A helpful contribution to this continuing effort would be an 
analysis of the heat flow and temperature gradient through composite walls 
for different operating conditions and composite wall proportions . An 
analysis of this type can establish the limits of usefulness of ceramic 
coatings. Accordingly, temperature -time histories were computed in a 
ceramic -coated metal-base composite wall by a numerical integration pro- 
cedure. Temperature -time histories in a ceramic-coated nozzle were ob- 
tained experimentally on a solid-propellant rocket in order to check the 
calculations . 

In the analytical study the effects of coating thickness, metal-base 
thickness, and chamber pressure on the maximum running time for a compos- 
ite wall were determined. The maximum running time was determined by 
limiting the local ceramic and metal temperatures to values lower than 
their melting points. 
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Experimental temperature -time histories of the metal section of the 
composite wall were determined for several chamber pressures and nozzle 
area ratios. A comparison of the agreement between the experimental and 
analytical temperature-time histories was made. 


APPARATUS AND PROCEDURE 


Motor 


A photograph of the solid-propellant rocket motor installation is 
shown in figure 1. The propellant was a composite, end -burning grain 
that was inhibited at the head end and circumferentially. An insulating 
sleeve was used to protect the chamber walls. The calculated gas equilib- 
rium chamber properties were as follows - 


Temperature , °R 

5122 

X 

1.164 

Cp, Btu/(lb)(°F) 

0-676 

c*, ft/sec 

4861 

Gas composition. 


mol percent 


Cl 

0.7 

H 

.3 

CO 

11.2 

C0 2 

11.6 

H 2 

5.5 

h 2 o 

39.4 

HC1 

20.8 

N 2 

9.0 

NO 

.1 

0 2 

.1 

OH 

1.5 


Nozzle 

A cross-sectional view of the nozzle geometry used in the experimen- 
tal investigation is shown in figure 2. The nozzle consisted of a copper 
insert at the throat and a steel divergent section. The nozzle throat 
area was varied by modifying the copper insert and the divergent section 
of the nozzle. Internal surfaces of the nozzle were sand-blasted, flame 
sprayed" with a thin (0.005 in. or less) coating of Nichrome, and then 
"flame sprayed" with zirconium oxide; the latter coating was sanded to 
make the surface and thickness more uniform. 
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Instrumentation 

Transient measurements of chamber pressure and. wall temperatures 
were recorded with a high-speed direct-recording oscillograph. Chamber- 
pressure static taps located just ahead of the nozzle measured essen- 
tially total chamber pressure. The method of installing wall thermocou- 
ples is illustrated in figure 2. Temperature measurements were made at 
nominal area ratios of 4.5 and 7.5. No temperature measurements were 
made in the thin ceramic-insulation coatings. Chamber-pressure and tem- 
perature measurements were accurate to ±1 percent. 

Procedure 

The weight of solid propellant burned was determined by weighing the 
motor before and after firing. Firings were made at several chamber 
pressures by varying the nozzle throat diameter. The measured metal wall- 
temperature histories were compared with those calculated for the same 
operating conditions . Symbols are defined in appendix A, and the method 
of calculating transient wall temperatures is presented in appendix B. 
Calculated values of the heat-transfer coefficient and the gas tempera- 
ture were determined as described in appendix C. The ceramic-coating and 
metal- wall thicknesses were measured before and after each min in the 
areas where instrumentation was located. There was no measurable erosion 
of the coating. The accuracy of the ceramic-coating measurements was 
limited to ±5 percent because of the inherent roughness of the coating. 


RESULTS AND DISCUSSION 

Factors Affecting Design of Ceramic-Insulated Nozzles 

An illustration of typical temperature histories calculated for a 
steel wall insulated with zirconium oxide at the nozzle throat is shown 
in figure 3. The insulation surface, interface, and metal cold surface 
temperatures are plotted against time for a 400-pound-per-square-inch- 
absolute chamber pressure. The rate of rise of the insulation surface 
temperature is very high at first and then abruptly decreases. The pri- 
mary reason for the decrease is that the high surface temperature reduces 
the heat transferred to the wall by decreasing the temperature difference 
between the gas and the wall . The quantity of heat that can be absorbed 
by the metal heat sink then determines the rate of insulation temperature 
rise . 


Optimization of insulat ion and metal thicknesses . - Since the heat 
absorbed by the nozzle wall is proportional to the temperature difference 
between the gas and nozzle wall, the higher the insulation surface tem- 
perature, the greater the reduction in absorbed heat. However, because 
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the melting temperature of ceramic insulating coatings is usually less 
than the flame temperature of high-performance propellants, the insula- 
tion temperature must be limited if severe erosion is to be avoided. 
Likewise, the metal base has a maximum temperature limit due to strength 
considerations. Therefore, in order to keep the nozzle weight to a mini- 
mum and stay within safe temperature limits on both insulation and metal, 
the insulation and metal thicknesses must be selected carefully for a 
given operating condition. 

An example of the effect of insulation and metal-base thickness on 
running time for the nozzle throat was calculated and is shown in figure 
4 for a chamber pressure of 400 pounds per square inch absolute and a gas 
temperature of 4660° F. The calculations were made for a composite wall 
consisting of zirconium oxide insulation (melting temp. 4500° F) on a 
mild-steel base. The limiting temperatures of the zirconium oxide and 
steel were taken to be 4100° and 1800° F, respectively. It can be seen 
that there is a well-defined optimum value of insulation- coating thick- 
ness that results in a maximum running time for each wall size. This is 
a result of the limiting temperatures of both the insulation and metal 
base. At the peak of each curve, both the insulation and the metal are 
locally at their limiting values. For larger than optimum insulation 
percentage, the insulation temperature reaches its limiting value first, 
while the metal temperature is limiting for insulation percentage less 
than optimum. It can also be seen that maximum nozzle life increases 
with increasing wall thickness or heat sink. However, the increase is 
not proportional to the wall thickness, because the larger temperature 
gradient that exists with heavier walls results in a lower average wall 
temperature, and the heavier wall is therefore a less efficient heat 
sink. The curves in figure 4 apply only to conditions at the nozzle 
throat. Primarily as a result of the variation in convective heat- 
transfer coefficient, the optimum required metal heat sink and insulation 
thickness will also vary through the nozzle. It is therefore necessary 
that the optimum proportion between insulation and metal thicknesses be 
determined at several points to minimize the wall weight along the entire 
nozzle. 

Effect of chamber pressure on maximum running time . - The rate of 
heat transferred to the nozzle walls directly affects the running time of 
a given nozzle. One of the important factors affecting heat-transfer 
rate over which the designer has some control is the chamber pressure. 

The convective heat-transfer coefficient and the chamber pressure are re- 
lated approximately as follows : 


The use of low chamber pressures is then a means of reducing the heat- 
transfer rate and increasing the running time of heat-sink nozzles. An 
example of the effect of chamber pressure on maximum nozzle running time 
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calculated for a nominal 0.4-inch composite wall, nozzle throat condi- 
tions, and a propellant gas temperature of 4660° F is shown in figure 5. 
The maximum nozzle running time was determined by optimizing the insula- 
tion thickness at each chamber pressure while holding the nozzle wall 
weight per unit area constant. This allowed the zirconium oxide insula- 
tion and steel base to reach their maximum limits of 4100° and 1800° F, 
respectively. It can be seen in figure 5 that, for low- chamber-pressure 
operation, very long running times are possible; while at high chamber 
pressures running time is severely limited. For example, a running time 
of 105 seconds was calculated for a chamber pressure of 50 pounds per 
square inch absolute , while at 1000 pounds per square inch absolute a 
running time of only 6.5 seconds is possible. While this trend results 
primarily from the reduced heat-transfer rate at lower chamber pressures, 
a second benefit is that a lower heat-transfer rate reduces the tempera- 
ture gradient in the nozzle wall, resulting in more effective use of the 
metal heat sink. It is also apparent that the optimum insulation and 
metal thicknesses are affected by variation in chamber pressure. The in- 
sulation layer comprises 45.83 percent of the wall thickness at 50 pounds 
per square inch absolute and only 4.50 percent at 1000 pounds per square 
inch absolute. 

While reductions in chamber pressure permit longer nozzle running 
times for a given wall weight per unit area, motor thrust is reduced. In 
order to keep motor total impulse constant (thrust multiplied by running 
time equals a constant), lower-chamber-pressure operation requires com- 
pensating increases in running time or nozzle throat area or a combina- 
tion of both. Increased nozzle throat area also results in a small reduc- 
tion in heat transferred to the nozzle due to reduced convective heat- 
transfer coefficients. It can be shown that the weight of a heat-sink 
nozzle for high -chamber -pres sure operation will be approximately the same 
as for a low-chamber-pressure nozzle with increased throat area and run- 
ning time to maintain constant total impulse . The advantage of low- 
chamber-pressure operation lies in the fact that a much greater range of 
running time is available for a required total impulse. It should be 
pointed out that low-chamber-pressure operation is best suited to high- 
altitude applications, since here the nozzle thrust coefficient is rela- 
tively unaffected by the resulting reduction in nozzle pressure ratio. 


Comparison of Calculated and Experimental Temperature Histories 

In order to determine how well the calculation procedure can predict 
temperatures in composite walls, experimental wall temperature measure- 
ments were made on a ceramic -insulated steel nozzle during solid- 
propellant rocket firings and were compared with calculated values. Ex- 
perimental temperature histories were obtained from fine thermocouples 
that were imbedded in the metal portion of the nozzle. The method of 
calculating average experimental heat-transfer coefficients and effective 
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combustion temperature Is described in appendix C. A comparison of a 
number of experimental and calculated temperature histories for the steel 
section of the composite wall at nominal area ratios of 4-5 and 7.5 and 
at several chamber pressures Is shown in figure 6. Examination of the 
temperature histories indicates that for all cases the measured metal tern 
peratures were lower than those calculated. The measured temperature rise 
for all the data varied from 71 to 86 percent (an average of 78 percent) 
of the calculated temperature rise at the end of the run. A primary rea- 
son that the ratio of the measured to calculated temperature rise varied 
from run to run was that the measurement of the insulating- coating thick- 
ness was only accurate to ±5 percent. 


One source of discrepancy between the measured and calculated tem- 
peratures was the assumption of one-dimensional (slab) heat flow. This 
is due to the fact that the annular nozzle wall has a larger cross- 
sectional area (or heat sink) than a flat wall of the same thickness. 

The effect of the assumption of slab heat flow on wall temperature can be 
estimated by assuming that local wall temperature is an inverse function 
of the wall cross-sectional area or, 

^slab ^annular 
^annular A slab 

and, if d 0 and &± n are the outside and inside diameter of an annular 
wall ; 

^annular 
A slab 

Applying this approximation to the data from figure 6, the difference be- 
tween slab and annular heat flow accounts for about 5 to 10 percent of 
the calculated temperatures. If the calculated temperatures are corrected 
to represent annular heat flow, the measured temperature then averages 
about 86 percent of the calculated temperature rise. 

Another factor that contributed to the difference between measured 
and calculated temperatures was the use of an end -burning solid-propellant 
grain to obtain the experimental temperature histories. The measured wall 
temperatures were probably decreased as a result of the local cooling ef- 
fect due to evaporation and ablation of the inhibitor and case insulation 
that encased the propellant grain. 



SUMMARY OF RESULTS 

1. Calculation of typical temperature histories in ceramic -insulated 
heat-sink nozzles indicated that there is a distinct optimum insulation 
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and metal heat-sink thickness that will minimize total wall thickness for 
a given operating condition and required running time. This optimum ex- 
ists when safe operating temperature limits are imposed on the insulation 
and metal. The optimum insulation and metal thicknesses will vary 
throughout the nozzle as a result of variation in heat-transfer rate. 

2. Reducing the chamber pressure results in a significant increase 
in the maximum running time of a given weight nozzle. Chamber pressure 
affects the optimum insulation and metal thicknesses . 

3 . Comparison of experimental and calculated wall temperature histo- 
ries showed that the measured temperatures averaged about 78 percent of 
the calculated temperatures for slab heat flow, and about 86 percent of 
the estimated temperature for annular heat flow. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 29, 1960 



APPENDIX A 


SYMBOLS 

area, sq in. 

characteristic velocity, ft/sec 

specific heat at constant pressure, Btu/(lb)(°F) 

diameter, in. 

ratio of nozzle outside to inside diameters 

convective heat-transfer coefficient, Btu/(hr)(sq ft)(°F) 

thermal conductivity, Btu/(hr)(sq ft)(°F/in.) 

dimensionless parameter, M = Ax^/(a At) 

dimensionless parameter, N = hf Ax/k 

number of slabs or increments 

Prandtl mmiber 

chamber pressure, lb/sq in. abs 
quantity of heat transferred, Btu/hr 
dimensionless parameter, R = Ax i k m /Ax m k i 
Reynolds number 
temperature, °F 

temperature at end of last wall increment, °F 
film temperature, °F 

adiabatic wall recovery temperature, °F 
surface temperature, °F 
flow rate, lb/sec 


thickness, in. 
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Ax incremental wall thickness, in. 

a thermal diffusivity, sq ft/hr 

X ratio of specific heats 

[i viscosity 

t time, sec 

At time increment, sec 

Subscripts : 
f film 

i insulation 

int interface between metal and insulation 

m metal 

s static 

t throat 

w wall 

1,2,3, .. . center of wall incremental thickness 
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APPENDIX B 


METHOD OF CALCULATING TEMPERATURE TRANSIENTS IN COMPOSITE WALLS 

A simplified numerical integration method was used to calculate 
temperature-time histories in a composite wall, using a high-speed digit- 
al computer. The method was based on the numerical method for calculat- 
ing one-dimensional transient heat conduction and storage described in 
reference 1. In this numerical integration procedure, the composite wall 
is divided into a number of increments Ax or slabs, as indicated in the 
following sketch: 


Composite wall cross section 


Insulation 


Metal base 



The temperature rise of a slab resulting from the heat absorbed during a 
finite period of time At is found by making a heat balance: 

Heat absorbed by slab = catering X Ax - q leaving X Ax 


where 


q =3 kA(dt/dx) 

and dt/dx ; the temperature gradient between adjacent slabs , is approxi- 
mated by assuming that the wall temperature between adjacent slabs varies 
linearly with x. 

Equations for transient conduction through a composite wall can be 
derived from the preceding heat balance. A detailed derivation of equa- 
tions will not be presented here but is available in reference 1 or 2. 
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The general form of the equations used to calculate temperatures at 
various points in the composite wall is as follows: 

I. The surface temperature t^, is given by 


21% + [M ± - (2N + 2)]t 0 + 2t-, 

t i - — — — i 

where the primed temperatures represent the temperature after an elapsed 
time At. The value of Mi was chosen so that Mi > ZN + 2 in order 
that tQ would not exert a negative influence on t^. The convective 
heat-transfer coefficient, which is contained in the parameter N, was 
found using the equation (ref. 2) h f = 0.023(k/d)Re°- 8 Pr 0 * 33 , where all 
fluid properties were evaluated at the film temperature [t f = (t s + t w )/2]. 

II- In the insulation layer 

Vi + ( M i - 2)t n + t n+ i 

*1,2,3, 4, 5 - ET 

III. At the interface between insulation and metal 

/ M i MmR \ 

( Unt-l + \2 + 2 "l - ^jUnt + ^int+l 

irrt ~ Mi jy^R 

~2 + ~ 2 ~ 


In this equation a choice of must be made, and there are several re- 

strictions governing this choice. The time increments must be equal 
(Arp = AT m ) in order to pass continuously in time increments from the in- 
sulation into the metal, and secondly the metal slab thickness A^ 
should be integral with the metal thickness x^. Therefore, in the equa- 
tion for Mm , 

_ (A *m )2 (*nAn ) 2 

% = ~ a m 


a number of slabs n^ was selected that made M m equal to or slightly 
greater than Mj_ . 

IV. In the metal section 


+ i 

D, 7, 


t n-l + (^m ~ 2 ) t n + Vt-l 
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V. The equation for calculating tg was 

Zt e-1 + (Mm ' 2 He 



The preceding calculation procedure has several assumptions or 
limitations : 

(1) One-dimensional heat flow (slab flow) 

(2) All heat absorbed in wall (perfectly insulated) 

(3) Constant average thermal diffusivity of insulation and metal 

( 4 ) Constant average gas film heat-transfer coefficient 

To begin the calculation procedure , an average surface temperature 
was assumed, and the convective heat-transfer coefficient was evaluated. 
An assumption of the average thermal diffusivity of the insulation and 
metal was then made. A choice of the number of insulation increments or 
slabs was made, and the dimensionless parameters N and M were then 

evaluated. The use of a large number of slabs improves the accuracy of 

the calculations , however, at the expense of an increase in the time re- 
quired to calculate the problem. The calculation procedure was repeated 

until the assumed surface temperature and thermal diffusivities agreed 

with the calculated values. 
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APPENDIX C 


METHOD OF CALCULATING EXPERIMENTAL VALUES OF 
LOCAL HEAT-TRANSFER COEFFICIENT 

One of the requirements necessary to calculate local heat-transfer 
coefficients is the knowledge of the effective flame temperature. The 
theoretical flame temperature was evaluated for the solid propellant used 
in the experimental investigation in the manner described in references 
3 and 4. The effective flame temperature was found by correcting the 
theoretical flame temperature for combustion inefficiency: 

* , | G measured 

^effective “ ^theor. I c * 

\ uheor. 


The recovery temperature was 


tR “ t g + °- 9 ( t effective " t s^ 

Local values of the gas film heat-transfer coefficient were calcu- 
lated using the following general relation: 


h f = 0.023(k/d)Re°- 8 Pr 0 ’ 333 


(ref. 2) where all properties were evaluated at the film temperature 
t^ = (t + t w )/2. By substituting and combining terms , the preceding 
equation can be rearranged into the following form: 


h 


f * 


0.023 


0.2 , , 

c p,fPf Im 

d°.2pr0. 67 \A t f 


The values of n, k, and Pr in this equation were estimated using the 
data from reference 5. The value of w/A, the rate of propellant flow 
per unit flow area, was determined in the following manner. The instan- 
taneous propellant flow rate for a solid-propellant motor can be approxi- 
mated very closely by 


^Pc 

/ Pc « 
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where w is the total propellant weight burned, p c is the chamber pres- 
sure , and j p c dt is the area under the chamber -pres sure trace . An 

average value of w for the time interval from ignition to tailoff was 
determined by substituting an average value of chamber pressure for the 
same time interval into the preceding equation. The value of heat- 
transfer coefficient calculated by this procedure is therefore an average 
value for the time interval from ignition up to but not including 
tailoff. 
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- Typical calculated temperature histories for composite vail at 
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Figure 4. - Calculated effect of insulation thickness and total wall thickness on nozzle 
throat running time. Maximum temperature limits of 4100° and 1800° F on ceramic coat in, 

and metal base. Chamber pressure, 400 pounds per square inch absolute; gas temperature 
4660 u F. 
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Figure 5. - Calculated effect of chamber pressure on maximum 
running time for nozzle throat for nominal 0.4- inch insulated 
vail. Gas temperature, 4660° Fj wall weight, 0. 1115 pound 
per square inch. 
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(a) Inner metal surface temperature. Measurement- station area ratio, a/A^, 4.76; 

insulation thickness, x i; 0.042 inch; base thickness, 3%, 0.355 inch; nozzle outside- 
to . inside-diameter ratio, c^/d^, 1.18; average chamher pressure, p c , 973 pounds per 
square inch absolute. 

Figure 6. - Comparison of typical calculated and measured temperature- time histories in 
steel nozzle vail insulated with zirconium oxide. 
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(c) Inner metal surface temperature. Measurement -station area ratio, a/a+, 4.43* 

insulation thickness, 0.033 inch; base thickness, x m , 0.353 inch; nozzle outside- 

to inside- diameter ratio, & 0 /&± n , 1.135; average chamber pressure, p p , 197 pounds per 
square inch absolute. c 

Figure 6. - Continued. Comparison of typical calculated and measured temperature- time 
histories in steel nozzle vail insulated with zirconium oxide. 
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Figure 6. - Continued. Comparison of typical calculated and measured temperature -time histories in steel nozzle 
wall insulated with zirconium oxide. 



Metal temperature 


25 



Running time, t, sec 

(f) Inner metal surface temperature. Measurement- station area ratio, A/A £, 7.51; insu- 
lation thickness, x-^, 0.030 inch; base thickness, x^, 0.359 inch; nozzle outside- to 
inside-diameter ratio, c^/d^, 1.097; average chamber pressure, p c , 197 pounds per 
square inch absolute. 

Figure 6. - Concluded. Comparison of typical calculated and measured temperature- time 
histories in steel nozzle wall insulated with zirconium oxide. 
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